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INTRODUCTION 


While  the  kinetic  theory  of  the  solid  end  gaseous  phases  of  materials 
are  fairly  well  developed,  little  is  known  about  the  kinetics  of  the  liquid 
phase.    The  study  of  the  velooity  of  the  interface  between  the  liquid  and 
solid  in  a  freezing  undercooled  liquid  was  undertaken  to  yield  information 
on  this  subject. 

When  a  sample  of  liquid  is  placed  in  a  long  tube  and  undercooled,  it 
is  observed  that  after  freezing  is  initiated,  the  freezing  front  moves  down 
the  tube  with  a  finite  and  reproducible  velocity.    The  velocity  is  directly 
dependent  on  the  degree  of  undercooling* 

In  this  experiment  an  attempt  is  made  to  develop  the  relationship 
betaeen  the  velocity  and  the  temperature  of  the  liquid.    The  results  were 
compared  to  a  similar  experiment  by  Barter  (l)  with  benzene  and  with  the 
theory.    An  empirioal  equation  is  proposed  whioh  fits  the  data. 

EQUIPMENT  AND  PROCEDURE 

The  sample  holder  was  a  closed  pyrex  glass  tube  with  a  oonoentric  well 
at  the  top  to  be  used  to  initiate  freezing.    Nitrobenzene  wee  vacuum  distilled 
into  the  sample  tube  through  a  side  arm  attaohed  to  the  top  of  the  sample  tube. 

The  sample  was  placed  in  a  oooling  bath  of  methanol  which  wee  held  in  a 
transparent  Dewar  Flask.    This  container  was  in  the  form  of  a  Dewar  Flask  to 
prevent  frosting  and  to  ensure  effioient  cooling.    The  outside  diameters  of 
the  tubes  making  up  the  Dewar  flask  were  51  mm  and  52  mm.    Plate  I  shows  the 
construction  of  the  sample  holder  and  bath. 

The  bath  was  cooled  by  passing  air  through  copper  coils  placed  in  liquid 
air  and  then  through  a  copper  U-tube  whioh  was  placed  in  the  oooling  bath. 


EXPLANATION  OF  PLATE  I 

Fig.  1    Cooling  Bath 
A*  Methanol 

Fig*  2    Sample  Tube 

B#    Nltro  benzene 
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PLATE  I 


The  temperature  of  the  sample  was  varied  by  adjusting  the  rate  of  the  flow 
of  air* 

The  temperature  of  the  bath  was  measured  using  a  copper- conatantan 
thermocouple  with  the  reference  junction  in  a  water-ice  mixture,  and  a  Rubicon 
potentiometer.    The  other  junction  was  sheathed  in  glass  and  plaoed  in  the 
cooling  bath.    A  stirring  apparatus  was  used  to  eliminate  temperature  dif- 
ferences in  the  bath.    The  pressure  was  maintained  at  720  ram  of  Hg  by  attach- 
ing the  side  arm  of  the  sample  tube  to  a  mercury  manometer  and  an  air  reservior. 

A  millimeter  scale  was  attached  to  the  outside  of  the  cooling  bath  parallel 
to  the  sample  tube.    The  position  of  the  freezing  interface  and  the  reading  of 
a  stopwatch  were  recorded  on  Kodak  Plus  X  Reversal  movie  film  using  a  Filmo  70 
Bell  and  Howell  movie  oamera  set  at  12.6  frames  per  second  and  placed  2.7  feet 
from  the  sample. 

Using  a  slide  projector  the  film  was  projected  on  a  screen  so  that  the 
image  appeared  about  1.5  times  as  large  as  the  actual  size  of  the  apparatus. 
Readings  were  taken  of  the  position  of  the  freezing  front.    Readings  of  tine 
were  taken  from  the  film  from  a  timing  devise  which  was  placed  near  the  sample 
and  soale  so  that  all  would  appear  in  the  picture. 

The  position  was  then  plotted  against  the  time  to  obtain  the  velocity. 
In  almost  all  oases,  the  points  fell  extremely  close  to  a  straight  line.  The 
slope  of  this  line  gave  the  velocity  in  each  particular  case.    The  data  are 
shown  in  Table  I. 


Turnbull  (8)  develops  a  theory  of  phase  boundary  migration  for  a  one 
component  system  where  the  effects  due  to  density  are  ignored*    oonsider  the 


DISCUSSION 


5 


Table  I,    Temperatures  in  Millivolts  and  Corresponding  Temperatures  in 
Degrees  Centigrade  and  Corresponding  Velosities  in  Centimeters 
Per  Second  * 


Temperature 

Temperature 

Velocity 

-  .289 

-  7,5 

6,3 

-  .169 

-  4,5 

fx  ry 

2,7 

-  .115 

-  3,0 

1,3 

-  .236 

-  6,3 

4,7 

-  .198 

-  5,2 

2,78 

-  .123 

—  3.2 

1,71 

-  .132 

-  3.5 

1,85 

-  .187 

«•  4.9 

3,10 

-  .272 

-  7.2 

6,2 

-  .070 

-  1,8 

0,94 

-  .206 

-  5.4 

3.7 

-  ,305 

-  8,1 

6,45 

-  .311 

-  8,2 

6,8 

-  .327 

-  8,7 

7,3 

-  .210 

-  5,5 

3,3 

-  .247 

-  6,5 

3,3 

-  .324 

-  8,5 

6,7 

-  .351 

-  9.2 

7,5 

-  ,050 

-  1.3 

0,7 

-  .039 

m  0,9 

0,59 

-  ,224 

-  5,9 

4,4 

-  .054 

-  1,5 

0,85 

-  ,252 

-  6,6 

5,8 

-  ,219 

-  5,8 

4,8 

-  ,259 

•  6,8 

5,35 

-  ,088 

-  2,5 

0,93 

-  ,278 

-  7,3 

6,8 

-  ,059 

-  1.6 

0.67 

-  .236 

•  6.2 

4,00 

-  ,151 

-  4.0 

2.45 

-  ,260 

-  7.1 

5.9 

-  .141 

-  S.7 

2.06 

-  .041 

-  1.0 

0,65 

-  8.7 

6,96 

-  .138 

-  3.3 

1.73 

-  .268 

-  7.0 

5.93 

-  .111 

-  2,9 

1,62 

*    Temperature  in  degrees  corresponding  to  that  in  millivolts  for  a  copper 
constantan  thermocouple  was  determined  from  data  by  Dike  (2), 
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oross  the  interface  from  the  liquid  to  the  solid  a  molecule  must  aquire  an 
activation  energy  AG-  •    The  frequency  of  transfer  is 

where  Vo  is  the  frequency  of  molecules  striking  the  barrier* 

If  ZXG-'is  the  difference  in  Gibbs  potentials  per  molecule  between  liquid 
and  solid,  the  frequency  of  transfer  in  the  opposite  direction  is 

Vsl  -  Vo  exp[-{AG-+AG-')/ft~[  J  . 
The  net  frequency  of  transfer  across  the  phase  boundary  is 

and  the  velocity  of  migration  of  this  boundary  is  proportional  to     V  • 
Then 

Frenkel  (4)  indicates  that  the  rate  of  change  in  phase  must  be  inversly, 
proportional  to  the  viscosity.    If  this  is  also  true  of  the  velocity  of 
crystallisation,  and  if        is  assigned  the  value  of  the  activation  energy 
associated  with  diffusion,  Turribull's  (8)  expression  becomes: 


where 


F(r)  -  l-expf-AG-'/feTl 
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If  be  written  as 


^1 


(where  G-<?in  order  that  &  O  at  T  =  To  )#  then  /~~  (  T)  oan  be  written 
as  a  Taylor' a  expansion  around  To 


where 


The  ooeffioients  of  F(t)  become i 


f"'/-r  >  =  fjZ  -  23l  +  aT3  +.  ±3!l  _  j  ^'a''/JtL 


Ik  To 

±5'  _ 

,    /*  3'3 
T    A*  To* 

/r1^ 

in 


r  (t&)  .~  J  jo  t0     *-z^'t/rX3  TTn 


re* 


I 


If    <3 '  is  the  only  important  derivative  of    ^  then 

F'"t%\  -  +-  kJl^Jk^-h  F'fr*)3 

-To1-  To 

(  t>  To*  5 

If  3"  is  the  only  important  derivative  of  3  then 

F  (To)    -  J?-—LT0_) 

1  o 

F(T„)  =  - 
' o 

If  the  only  important  derivative  of  ^  "then 

&u >  -  -  m^aki 

Tc 

and  so  on* 

If  ^iG-'is  written  at 


-  Aff  +  PAV, 
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■where  R  is  the  Helmoholtz  free  energy,  P  is  the  pressure  and  A  V  is  the 
ohange  in  volume  from  liquid  to  solid,  then  it  is  seen  that 

-  #  f%  £  )  +  FAV 

Where  -fL  and  fs  are  the  partition  funetions  for  the  respeotive  states. 
Then  if  P  Z\V  is  oonsidered  constant 

r/f  7" 


■f-  2/P 
and  so  forth. 

Eyring  (4),  treating  the  liquid  as  if  some  of  the  moleoules  had  solid- 
like degrees  of  freedom,  and  some  moleoules  had  gas-like  degrees  of  freedom, 
has  successfully  calculated  some  thermodynamic  properties  of  some  liquids. 

If  yk  is  the  volume  per  moleoule  in  the  solid  phase,    and  Vl  is  the  volume 

A/  7v 

per  moleoule  in  the  liquid  phase  then 

/V 

is  the  percentage  of  moleoules  of  the  liquid  whioh  have  gas-like  degrees  of 
freedom,  and 

A/  ~  A/(VL  -Vs)   =  ///J£ 
14  Vl 


is  the  percentage  of  molecules  of  the  liquid  whioh  have  solid  degrees  of 
freedom* 

Using  this  approach,  the  partition  function  for  a  substance  suoh  as 
benzene  or  nitrobenEene  is,  for  the  solid 

and  for  the  liquid  (including  free  rotation)! 

N  Vs  3  N 


Where  £$  is  the  energy  of  sublimation  at  the  melting  point,    &  is  the  Einstein 
characteristic  temperature  at  the  melting  point,  /Y\^  is  the  number  of  equilibrium 
sites  accessible  to  a  molecule  besides  its  most  stable  position.  -^2-  ia  the 
strain  energy  stored  in  the  system  as  the  molecule  shifts  to  one  of  the  /V7^ 
sites.    The  frequency  of  vibration  in  the  lattice  is    V    ,  and  /%SjCm  the 
rotator  oonstants  of  the  molecule,  for  example  >f  =     ^377- *£  XA  where  TA 
is  the  moment  of  inertia  around  the/f  axis. 
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Than 


sinoe  the  vibrational  terms  are  assumed  to  be  the  same*  By  straight  forward 
differentiation  with  respect  to  T  on  obtains i 

+/v*  f-  £At  (/-exp(-  %))-h^.f  (l-exp  c-j^r)) 
Vu       KZ  J 


II 


i  Vs_  £s  _  Vs  cfexp(^^r)  3  3fV<.-Vs)  ? 

Vs  2 R ret  (exP  (^rr)  +  ">h ( ^TWr))  +-cizexr( j4&Tr ) 

+  J—    4.  3  (Vl  -Vs)  7 
3-T z  Z  VLT^  J 
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Expressing  the  coefficients  of  Rr)  as  functions  of  the  coefficients  of^  gives 
,„     _  T&C3  _  £^    .  ,  ^2  _  <i  ^/  c*       /_£/  \?  1 

3C  C,x       IZC,*CX  _  17  C?  _  /  C,   \  v  "7 


Barter  (l)  shows  that /Vf? for  benaene  ia  of  the  form 

2  t  SI 

*  * 

The  following  are  values  of  constants  appearing  in  the  partition  function 
for  bensene  as  given  by  Eyring. 


*&  =  278.60K 

£i= 10,383.8  oal 

77,00  ml 
vi=  87,28  ml 

1,38047  X  10"16  erg/dog 
h~  6.624  X  10-27  erg  aeo 
hi  -6,0228  X  102s  molecules/mole 
8.31  X  107  erg/mole  deg 
>»  QC  r  2B3 
^=2.67 
4  =  .12 
Q-  80.24°K 
&-  55.2°K 
V  r  9  X  1012  CPS 
B  -=0.17165  cm"1 

Q  is  the  Debye  temperature  of  solid,  and  S  is  the  Debye  temperature  of 

the  liquid. 

These  values  give 

%'=C,  -  -  73.  S  ,  2  c2  =  x/o7 

This  yields 


F"(Tq)  r 
X 


-3  6,  &OOXJ0 


It  can  be  seen  from  the  form  of  the  expressions  for  fj  »  3  /  and  ?  " 
that  the  values  for  higher  order  derivatives  of  ^  beoome  progressively  less 
important  for  any  ;iven  set  of  constants*    Since  ^V-J^and  C,  are  negative* 
and  £;  is  by  far  the  most  dominant  term*   Eyring's  theory  will  give  a  seriea 

of  terns  of  alternating  sign*   Aa  the  seriea  1?  of  the  fora/MT -/^atz +A2at3  

and  A  is  equal  to  -  J  91  X  /o       »  it  oan  be  seen  that  this  form  does  not 
agree  at  all  with  empirical  results*    If  Ci  Is  positive  all  terns  are  nega- 
tive resulting  in  negative  velocities* 

RESULTS 

The  velocity  of  crystallisation  of  nitrobenseno  was  plotted  as  a  func- 
tion of  temperature  as  shown  in  Plate  II*   Barter  (1)  found  a  similar  plot 
for  benaene  as  shown  in  Plate  III*   Barter  obtained  an  empirical  equation 
for  benaene  in  the  form 

\Z-^(dATz  f  AT1  ) 
where  (if  at  =:  T-To ) 

and  yj  is  the  viscosity* 

By  extrapolating  viscosity-temperature  values  obtained  from  the  Handbook 
of  Chemistry  and  Physics  (6)  a  plot  was  obtained  of  VY\  versus  AT*  •  This 

is  shown  in  Plate  IV*   Using  these  values  and  obtaining  a  beat  fit  to  the 
original  data  an  empirical  equation  for  nitrobenzene  was  found  of  the  form 


V  -  A  (8  AtZ+AT'*) 


EXPLANATION  OF  PLATE  II 


Plot  of  velocity  of  orystallitatioa  of 
T'itrobergene  versus  temperature 


PLATE  II 


» 
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T  in  °c 


BHAttTin:  op  rum  in 


Plot  of  velocity  of  crystallisation  of 
Bensene  versus  temperature,  obtained  by  Barter 


PLATE  III 


EXPLANATION  OF  PLATE  IV 


A  plot  of  Vty  versus  ,£7"*  for  beneene. 
Values  of  V{  were  dote  mined  by  en  extrapolation 
of  values  from  the  Handbook  of  Chemistry  and  Physios* 


20 

PLATE  IV 
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where 

A  ~  .000  3  8? 

and 


This  is  shown  in  Plate  V.    As  oan  be  seen  by  comparing  the  plots  for 
benzene  and  nitrobenzene,  the  freezing  velocity  for  nitrobenzene  approaches 
zero  faster  as   T  approaches   To  •     This  leads  to  difficulty  in  measuring 
the  velocity  for  temperatures  close  to  7S  •    Hot  only  is  it  difficult  to 
start  freezing  and  keep  the  sample  at  an  equilibrium  temperature,  but  the 
freezing  front  is  poorly  defined  and  the  oontrast  is  not  great  on  the  film. 
In  some  instances  it  was  difficult  to  see  the  freezing  front  directly,  much 
less  on  the  film*    It  is  for  these  reasons  that  no  readings  takmi  corres- 
ponding to  supercooling  to  a  temperature  greater  than  zero. 

CONCLUSION 

The  combination  of  Turnbull' s  expression  with  Eyring's  partition  func- 
tion does  not  produce  agreement  with  the  empirical  results.    There  are 
several  places  from  which  this  discrepancy  might  arise. 

First,  the  expression  involving  the  relationship  between  the  frequency 
of  transfer  of  a  molecule  from  one  phase  to  another  to  the  activation  free 
energy  developed  by  Turnbull  might  be  incorrect.    Turnbull  assumes  there  is 
only  one  average  free  energy  in  each  state  or  only  one  barrier.    However  the 
activation  energy  could  depend  on  the  orientation  of  the  moleoules.    That  is. 


EXPLANATION  OP  PLATE  V 

Plot  of  velocity  of  crystallisation  of 
bensene  versus  temperature  using  expression 

V-  ^(BAT3-hA-T+) 


-   .  O  OO  387 


a 
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PLATE  V 


T  in  0  C 
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thj»re  may  b^i  a  rotation  barrier.    In  fact  there  might  be  a  great  number  of 

barriers  which  could  not  be  easily  represented  as  a  single  average  barrier. 

Ill  / 

1      Second,  in  taking  the  derivatives  of  AG  =AF+PAX/  with  respect  to  T 
PJ\  V  was  'assumed  to  be  constant,  or  at  least  the  change  in/^Vwas  assumed 
hegligable  when  compared  to  the  derivatives  of  A  F  .     This  assumption 
appears  to  be  well  justified  in  this  experiment. 

The  third  possible  source  of  the  discrepancy  between  the  theoretical  and 
empirical  results  is.fyrlng's  partition  function. 

Thug  last  appears  to  be  the  most  likely  souroe  of  difficulty  since 
Turnbull'''s  expression  seems  reasonable  and  the  higher  derivatives  of  A.  V/ 
diminish  much  faster  than  the  derivatives  of A  p  • 

Therefore /it  seems  likely  that  although  Eyring's  treatment  of  a  liquid 

j 

as  a  mixture-  of  moleojules  with  gas-like  degrees  of  freedom  and  moleoules 
with  solid-like  degrees  of  freedom,  while  accurate  enough  to  yield  many 
thermodynamics  properties,  is  not  sufficiently  accurate  to  give  correct 
values  for  the  velocity  of  crystallization. 
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A  sample  of  nitrobenzene  was  supercooled  to  various  temperatures  and 
freezing  initiated*    The  velocity  of  the  freezing  front  was  recorded  using  a 
movie  oamera  set  at  12.5  frames  per  second.    Readings  of  velocity  were  taken 
from  the  film  and  plotted  against  the  temperature.    These  results  were  com- 
pared with  a  similar  experiment  on  benzene* 

A  relationship  between  velooity  and  temperature  was  derived  using 
reaction  rate  theory  and  a  proposed  model  of  the  liquid  phase*    The  values 
obtained  did  not  agree  with  empirical  results*    An  empirical  equation  wbs 
obtained! 


